Relation between CPT Violation in Neutrino masses and mixings 
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Abstract 

The neutrino parameters determined from the solar neutrino data and the anti-neutrino pa- 
rameters determined from KamLAND reactor experiment are in good agreement with each other. 
However, the best fit points of the two sets differ from each other by about 10 -5 eV 2 in mass-square 
differenc and by about 2° in the mixing angle. Future solar neutrino and reactor anti-neutrino ex- 
periments are likely to reduce the uncertainties in these measurements. This, in turn, can lead to 
a signal for CPT violation in terms a non-zero difference between neutrino and anti-neutrino pa- 
rameters. In this paper, we propose a CPT violating mass matrix which can give rise to the above 
differences in both mass-squared difference and mixing angle and study the constraints imposed 
by the data on the parameters of the mass matrix. 
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I. INTRODUCTION 



KamLAND experiment has established the distortion due to oscillations in the anti- 
neutrino spectrum from reactors and determined the corresponding mass-square difference, 
A21, to a great precision 

a a. 

At present there is good agreement between A 2 i and 
the mass-squared difference of the neutrinos, A21, determined from the analysis of solar 
neutrino data 3, J]. However, the best-fit values of the two As differ from each other by 
about 10~ 5 eV 2 . Also, the best-fit mixing angles differ from each other by 2 to 3 degrees. 
Together, solar and KamLAND data impose the constraint |A 2 i — A21I < 1.1 x 10~ 4 eV 2 
js|. Future reactor experiments, located at a distance of about 70 Km from the source so 
that the oscillation minimum coincides with spectral maximum, are expected to improve 
the precision of anti- neutrino parameters even further 6|. Similarly future solar neutrino 
experiments [jj , are expected to improve the accuracy of neutrino parameters. If these future 
experiments confirm the present trend in the difference between neutrino and anti-neutrino 



aa. 



parameters, then CPT violation in the neutrino sector becomes an exciting possibility 

When it comes to the larger mass-squared difference, the atmospheric neutrino data 
prefers equal values for the neutrino and anti-neutrino mass-squared differences, though the 
uncertainties do allow large CPT violating effects [h]]. MINOS experiment is expected to 
measure the disappearance probability of muon anti-neutrinos with good precision in near 
future. If there is any difference between P[y^ — > v^) and Piv,,. — > v^), it will be a signal 



for CPT violation in the larger mass-squared difference also [111 ]. 

We assume that the main part of neutrino mass matrix is CPT conserving and it arises 
due to dynamics at an energy scale much below the scale at which CPT violation occurs. We 
further assume that CPT violation from a high scale, leads to an addition to the neutrino 
mass matrix of the form 

Mcpt = (1) 

where a and (3 are flavour indices and \i is a parameter with dimensions of mass. As we 
show in the next section, to reproduce the allowed differences between neutrino and anti- 
neutrino parameters, we require the scale of /1 to be ~ 10~ 6 eV. With this as input, we 
calculate the difference in the mass-squared differences and mixing angles for the neutrinos 
and anti-neutrinos. 

Quantum gravity effects can lead to CPT violation. The leading effective operators of 
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quantum gravity are suppressed as the inverse of the Planck mass Mp\. Such operators can 
give rise to CPT violating mass ~ v 2 /Mpi, where v is a low energy VEV of the quantum 
gravity model. If v = 174 GeV, the electroweak symmetry breaking scale, then we obtain 
fi ~ 10 -6 eV. The sign of the additional mass matrix for anti-neutrinos will have the opposite 
sign I2I ]. 

In eq. (CQ), A Q( g is a 3 x 3 matrix in flavour space. Quantum gravity effects are not sensitive 
to flavour. Hence it is expected that every term in the matrix X aj 3 is independent of both 
a and (3. We take this matrix to be of the form X a p = 1 for all a and f3. In this case, the 
CPT violating part of the neutrino mass matrix is of the form: 

//111. (2) 

V 1 1 1 / 

In our calculations, we take eq. (jSJ) as a perturbation to the main part of the neutrino mass 
matrix. The pattern of CPT violation in neutrino and anti-neutrino mass-squared difference 

n 

due to the above matrix was analyzed in [13J. Here we consider the constraints from data 
on the CPT violating parameters in mass-squared differences and mixing angles. 



II. CALCULATION 



We assume that the CPT conserving part of the light neutrino mass has real and non nega- 
tive eigenvalues Mj. In the mass eigenbasis, this matrix appears as M = diag(Mi, M 2 , M 3 ). 
We treat M as the unperturbed (0 th — order) mass matrix. Denoting the corresponding 
neutrino mixing matrix by U, we obtain the th — order mass matrix M in flavour space as 

M = U*MU f . 



(3) 



Explicitly, the matrix U has the form 



U 



( U el U e2 U e3 \ 

U^2 U^3 



(4) 



\ u Tl u r2 u t3 j 

where the nine elements are functions of three mixing angles and six phases. In terms of the 
above elements, the mixing angles are defined by 

U e2 



tan 9 



12- 



(5) 



u, 



fj3 



tan 023, (6) 
\U e3 \ = sm9 l3 . (7) 

In terms of the above mixing angles, the MNS matrix is written as 

U = diag(e if \ e i/a , e if3 )R(6 23 )A R(6 13 )A* R(6 l2 )dtag(e ia \e ia2 } 1). (8) 

The matrix A = diag(e^ , 1, e - ^) contains the Dirac phase 5. This phase leads to CP viola- 
tion in neutrino oscillations. a± and a 2 are the so called Majorana phases, which affect the 
neutrinoless double beta decay, fi, f 2 and f 3 are usually absorbed as a part of the definition 
of the charged lepton fields. It is possible to rotate these phases away, if the mass matrix in 
eq. ([3]) is the complete mass matrix. However, since we are going to add another contribu- 
tion to this mass matrix, these phases of the zeroth order mass matrix can have an impact 
on the complete mass matrix and thus must be retained. By the same token, the Majorana 
phases which are usually redundant for oscillations have a dynamical role to play now. 

Given the above Oth order (CPT conserving) neutrino mass matrix, we now add the CPT 
violating mass matrix to it. Thus the complete light neutrino mass matrix contains both 
CPT conserving and CPT violating terms. Given that fi is much smaller than the light 



neutrino mass scale (which should be greater than y/A atm eV), we can treat the CPT vio- 
lating mass matrix to be a perturbation of the CPT conserving mass matrix. The complete 
neutrino mass matrix in flavour space is 

M -> M' = M + fi\, (9) 

We assume that the symmetries inherent in M lead to tribimaximal mixing. But /xA breaks 
these symmetries. And hence the mixing angles given by the total mass matrix M' will not 
be tribimaximal. Below we compute the deviations from tribimaximality induced by /zA as 
well as the differences in mass-squared splittings. Note that these deviations will be equal 
and opposite for neutrinos and anti-neutrinos because /zA is CPT violating and is assumed 
to have opposite signs for particle and anti-particle. 

The perturbation formalism, by which the above computation can be done, was first 



developed in ref. 14j. Here we briefly recall the main features for completeness. The matrix 



relevant for oscillation physics is the following hermitian matrix 

M' f M' =(M + ^A) + (M + fiX). (10) 



To the first order in the small parameter fi, the above matrix is 

M t M + fiX^M + M^fiX. (11) 

This hermitian matrix is diagonalized by a new unitary matrix U . The corresponding 
diagonal matrix M' 2 , correct to first order in /i, is related to the above matrix by U' M' 2 U' . 
Rewriting M in the above expression in terms of the diagonal matrix M we get 

U'M' 2 U n = U(M 2 + m ] M + Mm)U\ (12) 

where 

m = fiU T XU. (13) 

Here M and M' are the diagonal matrices with neutrino masses correct to O^and l si order 
in /i. It is clear from eq. ffl~2l) that the new mixing matrix can be written as: 

U' = U(l+i5e), (14) 

where 50 is a hermitian matrix that occurs to first order in /i. Oscillation physics is un- 
changed under the transformation U — > UP, where P is a diagonal phase matrix. We can 
use this invariance to set the diagonal elements of the matrix 50 to be zero. 
From eq. ( fT2l) we obtain 



M 2 +m^M + Mm = M' 2 + [i5Q,M' 2 }. (15) 
Therefore to first order in /i, the mass squared difference AMfj = Mf — Mj get modified as 

Ail^ 2 = AM i + 2(M iJ Re[m ii ] - .U ; /.V,,, ., ). (16) 

The non-diagonal elements of 50 are given by 

iRejmij^Mj + Mj) /m(m tJ ) (M, - Mj) 

[ am; 2 am; 2 ' 1 j 



from which the changes in the mixing matrix can be computed by substituting 50 in eq. (1141) . 

The changes induced by the small parameter m are all proportional to the neutrino mass 
eigenvalues. They will have their largest values in the case of degenerate masses. Hence we 
assume degenerate neutrino masses Mj ~ M from hereon. In the expression for (50)^- in 
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Eq. (ITT]) , the second term is utterly negligible compared to the first, if we use degenerate 



15| 



masses. Thus we get a greatly simplified expression 

2iMRe(m ij 



(56), 



AMg 



(18) 



where we have substituted the Oth order mas-square difference in the denominator because 
the numerator already contains a factor of m. From Eq. ( ITS"]) , it is trivial to see that (50)\ 2l 
whose expression contains A 2 i in the denominator, is the largest among the (50)ij. 

Given the form of (50) ij, the elements of the modified mixing matrix can be obtained as 

U' aj = U aj + 5U a3 = U aj + i £ U ai (50) iy (19) 

i=i 

Knowing U' a j, we can define the modified mixing angles 9'^ in analology to the three equations 
given in Eqs. To compute 9'-, we first need to compute the changes in the five matrix 

elements 5U e j (j = 1,2,3) and 5U a $ (a = e, /i,r). Given that (50)i3, (£6)23 <C (50)i 2 , we 
can easily show that the changes in #13 and #23 are very small. To obtain 9' 12 , we need to 
evaluate the UM 



5U el 
5U e2 



-U e2 



Re(m 12 ) 
M 2 - Mi 
Re(m 12 ) 



(20) 
(21) 



M 2 - Mi 

For later convenience we define the complex numbers Zi = U e i + U^i + U T i, where U a i are, in 
general, functions of all six phases. 

In terms of the modified mixing matrix elements, 9 12 is defined as 



tan# 



12 



U, 



c2 



(22) 



Substituting the expressions from eqs. ([T9~j) - (j2T]) in eq. ff22|) . we get 

\xM I Z\ 1 1 z 2 j 



tan^ 



12 



tan 6»i2 + 2 



'AMfiCos 2 9 X2 



cos(ai + a 2 ) cos(<2i — a 2 ) 



tan 9u + e s 



(23) 



The modified solar mass-square difference is given by 



AM'i 



AM 21 + 2/iM \z 2 \ 2 cos(2a 2 ) - \zt\ 2 cos(2ai) 



AM 21 + e A 



(24) 
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Eqs. (1231) and (liMft give the modified mixing angle and mass-squared difference for neutrinos. 
The corresponding quantities for anti-neutrinos can simply be obtained by fi — > Thus 
we have 

AM' 2 l = A.Uj; • A A.U 2 - e A 
/ 

tan 6^ = tan #12 — eg = tan #12 — Eg (25) 

Note that the change in the mixing angle and the change in the mass-square difference 
have very different dependence on the Majorana phases a\ and a 2 . Therefore it will be 
straitforward to satisfy the experimental constraints for some combination of these two 
phases. 



III. RESULTS 



As mentioned in the introduction, we assume that the symmetries of the Oth order neu- 
trino mass matrix lead to tribimaximal mixing with #12 — 35.2°, #13 = and #23 = 45°. The 
)est fit value for the reactor anti-neutrino mixing angle, from KamLAND, is 9 l2 = 36.8° 
2|. For solar neutrinos it is 6' 12 = 32.6° [3], Thus we see that the anti-neutrino mix- 
ing angle is 1.6° more than the tribimaximal prediction whereas the neutrino mixing angle 
is 2.6° below the prediction. The differences between the best fits and the tribimaximal 
prediction are not equal and opposite. But, within the experimental uncertainties, they 
can be taken to be 2°. It is possible that the shifts in the mixing angles for neutrinos and 
anti-neutrinos are not equal. To explain such shifts, we need to invoke, along with CPT 
violating high scale physics, contributions from CPT conserving high scale physics, such as 
planck scale effects (l5j. The best fit value of reactor anti-neutrino mass-squared difference 



AM 21 = 8 x 10 5 eV 2 [2| and for solar neutrinos, the best fit value of the mass-squared 



difference is AM?' 2 = 6 x 10" 5 eV 2 



From eqs. (1231, (121ft and (125}. we find AMf 2 - AM|{ = 2e A = -2 x 10~ 5 eV 2 and 
6' 12 — 8 12 = 2ee = —4°. First we explore the following question: For what values of the CPT 
violating parameter /z, is there an agreement between the data and the hypothesis of CPT 
violating neutrino masses? In the introduction, we argued that fi ~ 10~ 6 eV. Here we take 
fi = p x 10~ 6 eV, where p is a number between 1 to 10, and derive the constraints the data 
imposes on p. We take degenerate neutrino masses for light neutrinos Mj = 2 eV, which is 
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the upper limit coming from tritium beta decay 17j and the neutrino mixing angles to be 
tribimaximal ones. Since #13 = in this case, the Dirac phase 5 can be set to zero without 
loss of generality. The zeroth order value of the smaller mass-square difference AM| X is 
set to 7 x 10~ 5 eV 2 , which is the average of the neutrino and anti-neutrino mass-squared 
difference. The phases /$ are set to zero. 

With these input values, the expressions for eq and e& become 



2e e = 0.04p[cos(2a 2 ) +cos(2ai)] = -4*tt/180 



2e A = 8p x 10" 



1 2 

- cos(2a 2 ) - - cos(2ai) 
o o 



2 x 10" 



Simplifying these equations, we get the following two conditions on cti and a 2 



(26) 



[cos(2a 2 ) + cos(2ai)] = — 1.75/p 
[cos(2a 2 ) — 2 cos(2oi)] = 7.5/p. 



(27) 



Solving these two equations and imposing the condition that —1 < cos(2ai), cos(2a 2 ) < 1, 
gives us the lower limit p > 3. For p — 4, eq. ( T2~7|) gives a\ = —70° and a 2 = 35°. For p = 6, 
these values change to a% = —60° and a 2 = 39°. Note that the Majorana phases a\ and a 2 
should necessarily be non-zero to satisfy the two constraints in eq. (|2"7|) . 



IV. CONCLUSIONS 



Both solar and reactor data are well explained by neutrino oscillations. Fit to solar data 
give a large region for the neutrino mass squared difference in the two flavor parameter space. 
The fit to reactor data however gives a very strongly constrained anti-neutrino mass squared 
difference. The best fits of the two mass squared differences are appreciably different from 
each other. Further improvement in KamLAND systematics and future solar neutrino data 
may further strengthen this discrepancy, thus giving a signal for CPT violation. We have 
demonstrated that flavour blind CPT violating neutrino masses from Planck scale physics 
can nicely accomodate this discrepancy, provided the Majorana phases of the neutrino mass 
matrix are appreciably large. This effect is crucially dependent on the neutrino mass spec- 
trum and gives rise to observable difference between A 2 i and A 21 only for a degenerate 
neutrino mass spectrum with m v ~ 2 eV, which is the largest allowed value from tritium 
beta decay data. The low value of the common mass implied by the WMAP bound [18( 
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leads to negligible difference between A21 and A21. This can however be compensated for 
by considering a slightly lower scale for the flavour blind CPT violating mass terms rather 
than the usual Planck scale. 

As we discussed in section II, the difference between A^ and A' 31 is negligible in this 
scenario if fi ~ 10~ 6 eV. If MINOS experiment were to observe a signal for CPT violation 



the above difference should be of order 10 3 eV 2 . Accounting for such a large CPT 



violation in the current scenario requires the CPT violating mass parameter to be of the 
order of 10 -3 eV. To obtain such a large value, the scale of CPT violating physics has to 
be three orders below the Planck scale. The flavour matrix X a p, can not be flavour blind 
because it would lead to an unacceptably large CPT violation for A21. Hence an appropriate 
texture should be imposed on \ a p. 
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